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过 冷水 滴 凝 固 特性 的 计算 与 实验 研究 


HCH, HER, Stim, Et, Wit, ER 

(中 国 空气 动力 研究 与 发 展 中 心 空气 动力 学 国家 重点 实验 室 ， 绵 阳 621000) 
摘 E: 针对 飞机 结 冰 的 特殊 凝固 过 程 ， 基 于 Enthalpy-porosity 模 型 ， 发 展 了 改进 的 结 冰 特性 预测 模型 及 数值 计算 
方法 ， 并 基于 凝固 实验 开展 了 所 建 模型 及 预测 方法 的 验证 。 研 究 表明 ， 相 对 于 传统 Enthalpy-porosity 方 法 ， 所 发 展 的 
方法 可 有 效 表 征 非 平 衡 凝固 过 程 ， 因 而 对 冷水 滴 凝 固 特性 的 表征 具有 较 好 的 改进 。 基 于 所 建 方法 ， 开 展 了 水 滴 凝 四 
特性 的 影响 因素 分 析 ， 获 得 了 不 同 过 冷 条 件 下 水 滴 凝 固 过 程 的 温度 分 布 及 相 界 面 变化 特征 。 相 关 研 究 可 为 结 冰 热力 
学 模型 的 改进 ， 以 及 结 冰 特性 的 精细 化 预测 提供 重要 参考 。 
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中 图 分 类 号 : V211. 3 文献 标识 码 : A 


Computational and Experimental Study on Freezing Characteristics of 


Supercooled Water Droplet 


XIAO Guang-Ming DU Yan-Xiao Li Wei-Bin WANG Qiao YI Xian GUI Ye-Wei 
(State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China) 


Abstract: Based on the Enthalpy-porosity model, an improved icing prediction model and numerical 
method are developed for the special freezing process of aircraft icing. The model and prediction method 
are validated by the related experiments. Compared with the traditional Enthalpy-porosity method, the 
results show that the new method can characterize the non-equilibrium freezing process of supercooled 
water droplet more effectively. Using the developed method, the influencing factors of the freezing 
characteristics of supercooled water droplet are analyzed, the temperature distribution and phase interface 
characteristics are obtained under different cooling conditions..The results can provide an important 
reference for the improvement of icing thermodynamic model and finer prediction of icing accretion. 
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引言 的 影响 特征 。 
近年 来 ， 随 着 结 冰 预 测 精度 要 求 的 提高 ， 结 冰 
飞机 结 冰 是 过 冷水 滴 撞 击 飞 机 机 体 表 面 并 发 生 。 ”过 程 的 过 冷水 滴 的 非 平衡 凝固 效应 引起 了 研究 者 越 
固 的 特殊 液 / 固 相 变现 象 ,也 是 影响 飞行 安全 的 重 。 来 越 多 的 关注 oa。 但 由 于 凝固 过 程 的 复杂 性 ， 目 
隐患 之 一 上。 在 飞机 结 冰 过 程 中 ,过 冷水 滴 的 凝 ”前 非 平衡 凝固 特性 预测 手段 的 研究 还 较为 缺乏 ， 巴 
通常 分 为 两 个 典型 阶段 外。 第 一 阶段 由 形 核 开 始 ，“。 测 精度 也 有 待 于 提高 ， 因 而 过 冷水 滴 凝 固 规律 特征 
水 滴 从 热力 学 非 平衡 态 过 渡 到 热力 学 平衡 态 的 阶 ”的 相关 研究 也 较为 薄弱 。 本 文正 是 针对 飞机 结 冰 过 
， 即 梳 蝇 形成 阶段 ， 也 有 研究 者 称 其 为 部 分 凝固 。 ” 程 过 冷水 滴 凝 固 的 非 平衡 效应 ， 发 展 了 过 冷水 滴 北 
段 四 。 第 二 阶段 为 完全 凝固 阶段 ， 即 液 固 相 界 固 特性 预测 方法 ， 并 自行 搭建 了 实验 系统 开展 了 所 
液 相 相向 固 相 推 进 直 至 凝固 完成 的 过 程 。 在 飞 ” ” 建 方 法 的 实验 验证 。 相 关 研 究 可 为 结 冰 热力 学 模型 
结 冰 过 程 中 ， 由 于 撞击 引起 的 异 相形 核 过程 ， 使 《的 改进 ， 以 及 结 冰 特 性 的 精细 化 预测 提供 参考 。 


固 的 第 一 阶段 显著 快 于 第 二 阶段 后 。 由 于 结 冰 物 


1 非 平衡 凝固 模型 


过 程 的 复杂 性 ， 因 而 目前 大 多 关于 过 冷水 滴 结 冰 
研究 把 重点 放 在 凝固 第 二 阶段 ， 即 相 界 面 的 推进 实验 研究 表明 ， 过 冷水 滴 的 凝固 过 程 由 两 个 典 


过 程 上 5。 鉴于 此 ,大 多 研究 者 基于 平衡 凝固 的 相 型 阶段 构成 。 在 第 一 阶段 即 枝 品 形成 阶段 ， 过 冷水 


关 
完 


理论 与 方法 开展 结 冰 特 性 的 预测 研究 ， 如 不 少 有 滴 温度 由 过 冷 态 上 升 至 凝固 点 ， 并 伴随 潜 热 的 部 分 


者 [&.9 基 于 Enthalpy-porosity 法 研究 凝固 过 程 的 液 释放 ， 是 非 平 衡 凝 固 阶 段 ， 第 二 阶段 为 由 热 扩散 驱 


相 变 行为 ， 但 该 方法 主要 针对 平衡 凝固 过 程 ， 无 动 的 界面 推进 阶段 ， 也 是 在 等 温 下 进行 的 平衡 凝 
法 表征 过 冷水 滴 凝 固 的 第 一 阶段 即 非 平衡 凝固 过 程 过 程 。 尽 管 研究 表明 凝固 第 一 阶段 相对 于 第 二 阶段 
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的 时 间 相 对 较 短 ， 但 由 于 结 冰 条 件 的 不 同 ， 使 第 一 
阶段 结束 时 形成 了 第 二 阶段 的 不 同 初始 条 件 ， 从 而 
影响 了 凝固 的 后 续 特 征 。 因 此 ， 过 冷水 滴 整 个 凝固 
过 程 的 有 效 预 测 应 综合 考虑 第 一 阶段 的 非 平 衡 凝固 
效应 及 第 二 阶段 的 平衡 凝固 效应 。 针 对 过 冷水 滴 的 
凝固 两 个 典型 阶段 的 特点 ， 如 何在 相 界 面 推进 过 程 
的 预测 中 考虑 第 一 阶段 的 非 平衡 效应 ， 是 有 效 预 测 
结 冰 全 过 程 凝固 特性 需 解决 的 重要 问题 。 


图 1 过 冷水 滴 凝 固 的 两 个 典型 阶段 


Fig.1 Two typical freezing stages of supercooled water 


droplet 
考虑 到 凝固 第 一 阶段 完成 并 形成 混合 态 的 过 程 
中 ， 过 冷水 滴 已 经 有 了 潜 热 的 部 分 释放 ， 因 此 ， 凝 
固 第 二 阶段 预测 中 应 进行 考虑 凝固 第 一 阶段 特性 参 
数 的 修正 。 将 无 量 纲 过 冷 度 & 表示 为 : 


E=cAT/L a) 


其 中 ，c, 为 固 相 比热容 ， 世 为 固 / 液 相 变 潜 热 ，AT 
的 计算 表达 式 为 : 


AT = T, = de Q ) 
其 中 ，7 为 液 相 凝固 温度 ， 克 ,ws 为 过 冷 态 温 


度 。 
凝固 第 一 阶段 结束 时 混合 态 中 未 发 生 相 变 的 液 


相 分 数 上 可 表示 为 中 


fabia (3) 
C 


Ss 


其 中 ，c 为 固 相 比热容 。 
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因此 ， 混 合 态 的 液 / 回 相 变 潜 热 Le 可 表示 为 : 


Lin =L f (4) 


混合 态 的 有 关 物 性 参数 PM. 可 表示 为 : 


PM,,,, = PM, f, +PM,(- f.) (5) 


其 中 ，PM, 和 PM , 分别 表示 对 应 的 液 相 和 固 相 物 


性 参数 。 

当 获 得 凝固 第 一 阶段 凝固 过 程 的 物性 参数 后 ， 
考虑 过 冷 阶段 非 平衡 凝固 效应 的 相 变 问题 即 可 转化 
为 平衡 凝固 问题 。 为 借鉴 平衡 凝固 条 件 下 液 / 固 相 变 
传 热 的 预测 方法 ， 作 以 下 几 点 假设 : 

(1) 忽略 水 滴 凝 固 过 程 水 滴 的 变形 ; 

(2) 由 于 撞击 形成 的 异 相形 核 作用 ， 凝 固 第 一 
阶段 时 间 显著 小 于 第 二 阶段 ; 

(3) 凝固 过 程 的 液 相 区 实质 为 已 在 凝固 第 一 阶 
段 发 生 了 部 分 相 变 的 混合 态 ， 因 此 凝固 第 二 阶段 计 
算 过 程 液 相 区 物性 参数 由 混合 态 参数 代 蔡 。 

基于 上 述 假设 并 借鉴 Enthalpy-porosity 方法 ， 
则 凝固 过 程 的 流动 及 传 热 控 制 方程 可 描述 为 : 


2 
Ox, 


Ou, Ou, Op O | Ou, 
p| —+u, = 一 + 
ôt Ox, Ox © Ox, | Ox; (7) 


—p8,B(T-T)+S, 
pfd, u T) 2 (a0 7 
ôt ’ Ox, ôx, | Ox, 


SE, mo u Aixo x 分 别 表示 不 同方 向 的 速度 


和 坐标 《对 于 三 维 问题 ， j=1,2,3)。p 和 了 分 


Oo 


别 对 应 于 压力 和 温度 变量 , T) 为 满足 Boussinesq 近 


似 下 的 参考 温度 ，h KAR. pr ur PADI 
对 应 于 密度 、 粘 性 系数 、 热 膨胀 系数 和 热 导 率 等 物 
性 参数 。g, 为 不 同方 向 的 重力 加 速度 ，$; 为 凝固 第 


二 阶段 两 相 区 的 源 项 ， 可 描述 为 中 : 
S,=-C(1-f,)u, (9) 


其 中 ，C 为 表征 凝固 过 程 两 相 区 特征 的 参数 ， 为 
结 冰 过 程 凝固 第 二 阶段 液 相 率 的 比例 ， 可 表示 为 : 


0 Per, 
T-T 
h= < TF aT ST (10) 
T-T, l 
1 T >T, 


其 中 ,和 工分 别 为 液 / 固 相 变 过 程 完全 熔化 和 完全 


凝固 时 的 温度 。 

采用 有 限 体积 法 离散 控制 方程 ， 非 稳 态 项 采用 
时 间 一 阶 精度 格式 ， 采 用 强 隐 式 迭代 法 进行 求解 。 
基于 固定 网 格 技术 ， 不 直接 追踪 两 相 界 面 ， 通 过 温 
SMW KA, REAANZI RM, ae 
免 直接 追踪 两 相 界面 带 来 的 困难 。 基 于 上 述 方法 ， 
对 项 目 组 前 期 自行 编制 的 相 变 传 热 计算 代码 鸣 进 
行 了 方法 的 改进 并 开展 了 典型 工 况 的 算 例 分 析 。 
2 非 平衡 凝固 模型 

以 直径 为 4mm、 高 度 为 2.5mm 的 过 冷水 滴 为 
例 ， 分 析 其 凝固 过 程 的 温度 分 布 及 相 区 分 布 特性 。 
其 中 , 水 滴 底 面 为 冷却 面 。 计算 中 底部 为 等 温 边 界 ， 
表面 采用 对 流 边界 。 计 算 物 性 参数 如 表 1 所 示 。 

表 1 ， 冰 和 水 的 材料 物性 参数 5 
Table 1 Physical properties of water and ice! 


HIE p ikem? clkgtK! A /Wem'K! L/k]-kg! 


DT 


J 


水 1000 4187 0.6 334624 
冰 920 2031 2.3 


图 2 显示 了 过 冷水 滴 凝 固 过 程 的 温度 分 布 特 
性 ;图 3 显示 了 过 冷水 滴 凝 固 过 程 的 相 区 分 布 特性 。 
在 凝固 初始 阶段 ， 相 界面 与 等 温 线 与 底部 冷 表 面 近 
似 平行 。 随 着 凝固 过 程 的 推进 ， 相 界面 与 等 温 线 逐 
步 向 水 滴 顶 部 弯曲 。 
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2 考虑 过 冷 效应 的 水 滴 凝 固 过 程 温 度 分 布 特性 


Fig.2 Temperature distributions of droplet freezing 


considering supercooling effect 
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3 ”考虑 过 冷 效应 的 水 滴 凝 固 过 程 相 区 分 布 特性 


Fig.3 Phase interface distributions of droplet freezing 


considering supercooling effect 
图 4 显示 了 过 冷 度 为 -10'C 条 件 下 考虑 与 不 考 
虑 非 平衡 凝固 效应 的 凝固 速率 比较 ; 图 5 显示 了 水 
滴 不 同 过 冷 度 对 凝固 相 界 面 变化 速率 的 影响 特征 。 
可 以 看 出 ， 在 考虑 过 冷 期 间 非 平衡 凝固 效应 影响 条 
件 下 ， 凝 固 速 率 较 传统 方法 有 所 提高 。 这 是 由 于 所 
建 模型 考虑 了 过 冷 阶段 形 核 和 村 晶 生长 过 程 中 水 滴 
潜 热 的 部 分 释放 现象 ， 因 而 在 界面 推进 阶段 ， 液 / 
固 相 变 潜 热 的 释放 量 相对 减少 ， 因 而 界面 推进 速度 
较 传 统 方法 有 所 提高 。 从 图 5 可 以 看 出 ， 在 相同 条 
件 下 , 水 滴 过 冷 度 越 大 , 相 变 的 驱动 力也 相对 越 大 ， 


NA y < _—_ 
因而 凝固 速率 也 相对 越 高 。 
1.2 
1 
一 » 
E 0.8 n" 
8 > a = 
£ 0.6 py 2 
2£ = 
E > ž m 
| w x 
2 a 
a . 
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图 4 过 冷 效 应 对 相 界 面 变化 速率 的 影响 


Fig.4 Effect characteristics of supercooled effect on freezing 


rate of phase interface 


Phase Interface (mm) 


0 15 20 


10 
Time (s) 
图 5 过 冷 度 对 凝固 相 界 面 变化 速率 的 影响 


Fig.5 Effect characteristics of supercooling on freezing rate 


of phase interface 

图 6 显示 了 水 滴 凝 固 过 程 不 同位 置 的 温度 变化 
曲线 。T01-T05 为 从 底面 到 表面 位 置 点 。 可 以 看 出 ， 
越 靠近 冷却 面 ， 温 度 降低 速率 越 高 。 这 是 由 于 凝固 
第 二 阶段 主要 是 在 相 变 温度 下 由 热 扩 散 驱 动 而 发 生 
等温 平衡 凝固 过 程 ， 导 热 是 该 过 程 的 主导 传 热机 
制 。 因 此 ， 越 靠近 冷却 面 ， 导 热 热 阻 相对 越 小 ， 温 
降 速率 也 相对 越 大 。 


Temprature (°C) 


0 20 40 


图 6 水 滴 不 同 高 度 位 置 的 温度 变化 特性 


Fig.6 Temperature variations of droplet freezing at different 


heights 

3 非 平衡 凝固 模型 
3.1 ”实验 系统 设计 

为 了 验证 所 发 展 模型 及 方法 的 有 效 性 ， 本 项 目 
自行 搭建 了 实验 系统 并 开展 了 相应 的 实验 研究 。 实 
验 系统 示意 图 如 图 7 所 示 ， 由 半导体 制冷 系统 、 
Agilent 34970A 多 点 数据 采集 仪 、ANV TF100 PID 
温度 监视 器 、 温 度 控制 器 、MotionXtra HG-100K 高 
速 摄像 机 、FLIR E60 红外 相机 、LED 无 影 光源 系统 、 
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工控 机 及 调 压 电源 组 成 。 实 验 时 由 半导体 制冷 台 提 
供 结 冰冷 环境 ， 由 PID 温度 控制 系统 将 制冷 台 表 面 
温度 维持 在 实验 所 需 的 温度 条 件 ， 基 于 界面 追踪 的 
方法 ， 采 用 高 速 摄像 机 记录 获得 水 滴 凝 固 过 程 相 界 
面 随时 间 的 变化 特性 ; 采用 红外 相机 记录 水 滴 凝 固 
过 程 温 度 的 实时 变化 特性 。 


Led Shadowless .©@ 


on 
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图 7 水 滴 凝 固 特性 实验 系统 示意 图 

Fig.7 Experimental setup for droplet freezing 
3.2 ”实验 凝固 特性 与 计算 结果 的 比较 

在 水 滴 凝固 实验 中 可 以 发 现 ， 在 过 冷 期 间 的 形 
核 与 枝 唱 生长 阶段 ， 水 滴 在 过 冷 条 件 的 驱动 下 逐步 
由 透明 态 向 模糊 态 过 渡 ; 当 水 滴 达 到 平衡 温度 ， 固 / 
液 界面 开始 出 现 并 由 液 相 区 向 固 相 区 移动 。 如 果 以 
固 液 /界面 出 现 的 起 始 时 刻 作 为 凝固 第 二 阶段 的 起 
点 ， 采 用 高 速 相 机 记录 相 界 面 随时 间 的 变化 特征 ， 
并 同时 采用 红外 热 像 仪 记录 水 滴 不 同时 刻 的 温度 分 
布 规律 。 图 8 显示 了 过 冷 度 为 -10C 和 条件 下 凝固 第 二 
阶段 相 界 面 随时 间 的 变化 过 程 ， 图 9 显示 了 相对 条 
件 下 相 界 面 推进 过 程 水 滴 红外 温度 场 随 时 间 的 变化 
特征 。 可 以 看 出 ， 在 凝固 第 二 阶段 ， 在 热 扩 散 的 作 
用 下 ， 相 界面 由 液 相 区 向 固 相 区 推进 ， 而 水 滴 的 低 
温 区 也 逐步 由 底部 向 顶部 发 展 。 
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图 8 水 滴 凝 固 过 程 相 界 面 随时 间 的 变化 特征 


Fig.8 Experimental evolution characteristics of phase 


interface 
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图 9 ”水滴 凝 固 过 程 温度 随时 间 的 变化 特征 
Fig.9 Experimental temperature variations of droplet freezing 
图 10 显示 了 的 计算 与 实验 相 界 面 移动 特性 的 
比较 ; 图 11 则 显示 了 过 冷 度 分 别 为 -10'C 和 -14'C 条 
件 下 ， 水 滴 表 面 计算 温度 与 实验 红外 温度 变化 特性 
的 比较 。 从 图 10 可 以 看 出 ， 在 相同 过 冷 度 条 件 下 ， 
改进 预测 模型 的 计算 结果 与 实验 结果 的 吻合 程度 优 
于 传统 Enthalpy-porosity 预测 模型 。 比 较 而 言 ， 考 
虑 过 冷 效应 的 改进 模型 预测 凝固 速率 要 高 于 传统 模 
型 。 从 计算 与 实验 结果 的 比较 可 以 看 出 ， 计 算 与 实 
验 相 界面 随时 间 的 移动 特性 吻合 较 好 ， 计 算 与 实验 
温度 变化 历程 虽 存 在 一 定 偏差 但 变化 趋势 基本 相 
似 ， 表 明 所 建 模型 的 有 效 性 。 
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Fig. 10 Comparison of experimental and simulated freezing 


rates of phase interface 
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图 11 水 滴 凝 固 过 程 实验 与 计算 温度 变化 特性 的 比较 


Fig. 11 Comparison of experimental and simulated 


temperature variations of droplet freezing 


4 结论 


本 文 针 对 飞机 结 冰 过 程 过 冷水 滴 结 冰 的 非 平衡 
凝固 效应 ， 开 展 了 过 冷水 滴 结 冰 凝 固 特性 的 计算 与 
实验 研究 。 主 要 结论 如 下 : 

(1) 基于 Enthalpy-porosity 模型 ， 发 展 了 过 冷 
水 滴 结 冰 特 性 预测 模型 及 数值 计算 方法 ， 并 自行 搭 
建 了 实验 台 ， 开 展 了 数值 方法 的 实验 验证 。 相 对 于 
传统 方法 ， 所 发 展 的 过 冷水 滴 结 冰 特 性 预测 模型 及 
数值 方法 能 有 效 表征 非 平 衡 凝 固 效 应 ， 因 而 可 用 于 
过 冷水 滴 凝 固 特 性 的 预测 ; 

(2) 在 结 冰 过 程 中 ,过 冷 度 越 大 , 水滴 的 凝固 
速率 相对 越 高 .但 在 考虑 非 平衡 凝固 效应 的 条 件 下 ， 
过 冷水 滴 凝 固 速 率 要 高 于 不 考虑 非 平 衡 凝 固 效应 的 
工 况 。 因 此 ， 飞 机 结 冰 过 程 中 过 冷水 滴 的 凝固 特性 
具有 与 平衡 凝固 过 程 不 同 的 特征 。 有 效 表征 结 冰 过 
程 的 非 平 衡 凝固 现象 是 改进 结 冰 热 力学 模型 ， 提 高 
结 冰 特 性 的 重要 基础 。 
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